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Abstract 

We have discussed conditions under which probabihty of oscihation (fg — )• f^) is 
independent of CP violating phase 5. The condition of magic basehne on its length 
is well-known. We have proposed another condition which is on neutrino energy. We 
have shown that magic baseline condition is not possible in general, for small ^13 
with non-standard interaction and for large ^13 with both standard and non-standard 
interactions. However, neutrino energy condition is possible for such cases as well as 
for cases where magic baseline condition is applicable. We have discussed how one 
may resolve hierarchy problem for neutrino masses by using such energy condition. 
For a baseline of length 650 Km, using this energy condition we discuss the possible 
number of fj,~ events at the detector for a period of 5 years and also the sensitivity in 
measurement of cos^ ^13. 
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The probability of oscillation of different flavors of neutrinos depends on various parame- 
ters present in the neutrino mixing matrix - the PMNS matrix [1] as well as the mass squared 
differences. Although two angles $12 and 623 are known with some accuracy but there is only 
upper bound for ^13 [2] and CP violating phase 5 is unknown. Although mass squared differ- 
ences for different neutrinos are known but the sign of Amsi^ (where Amf^- — m^ — rrij and 
is the mass of i-th neutrino) is unknown. Due to correlations among these unknowns there 
are ambiguities [3] in analysing neutrino oscillation data. To reduce such ambiguities it is 
useful to choose suitable baseline [4]. Particularly, magic baseline [5, 6, 7, 8] is useful for some 
specific length for which the perturbative expression of probability P^^^u^ is independent of 
S upto order (where a = Am^i/ Aml^). Then it is easier to find out other parameters 
apart from 6 which plays the role in neutrino oscillation. However, using magic baseline has 
its limitations also - namely (a) it may not be always possible to place the detector at a 
magic baseline distance from the source of neutrino production, (b) if there is non-standard 
interaction (NSl) [9] then we have shown that it is difficult to get 6 independent probability 
Pu^^Vf, using the magic baseline condition. 

In this work we present another condition - which may be termed as magic energy con- 
dition under which also the probability of oscillation will be independent of S upto order a^. 
Considering this condition one might be able to circumvent the above two shortcomings of 
the magic baseline. Unlike magic baseline condition this condition depends on length of the 
baseline also apart from its dependence on ^J^Gpfie (Here rif, is the electron number density 
of the matter). To use the magic energy condition, one option could be to analyse oscillation 
data in small energy bins. However, the other better option could be to use monoenergetic 
neutrino beam at the source as proposed in recent years [10, 11]. The idea is about using 
nucleus which absorbs an electron and emits a neutrino. By acelerating the mother nuclei 
with suitable Lorentz boost factor one may get the suitable neutrino energy which is satisfied 
by magic energy condition. Due to the monoenergetic nature of the beam it is expected to 
have better precision in finding various neutrino oscillation parameters. 

The magic baseline condition was initially obtained using the perturbative expansion for 
small Qxz with Standard Model interaction (SMI) [5]. Here, we obtain the modified form of 
magic baseline condition for both small and large ^^13 and for both SMI and NSI. Besides, we 
obtain the magic neutrino energy condition for both SMI and NSl and also for both small 
and large ^13 as allowed by present experiment [2]. Finally, we have compared advantages 
and disadvantages in considering magic baseline condition and magic energy condition in 
experiments. 
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Flavor eigenstates may be related to mass eigenstates of neutrinos Vi as 

>= ^^/mli^i >, U = R22.Ri2.{5)Ri2 and i = 1,2,3, (1) 

i 

where U is PMNS matrix [1] and Rij are the rotation matrices. General probabihty expres- 
sion for oscillation of neutrino of flavor / to neutrino flavor m in matter (satisfying adiabatic 
condition for the density of matter) is given by 

P{n ^ i^m) = - 4 MJ'iT] Sin' A;^. + 2 Im[f^] sin 2 A;^. (2) 

i>j i>j 

where J^J" = U'l^U'i-U'^^U'^j and A-^- = A'm|L/(4£;). Here A'mfj = -m^ and label ( ' ) 
indicates the neutrino matter interaction induced quantities corresponding to those quantities 
in vacuum. Let us write x = A'^^ , y — A32 and z = A'i2- Using trigonometric identities : 
— sin^ x+sin^ y— sin^ z = 2 sin x cos y sin z and — sin 2a;-|-sin 2y — sin 2z — —4 sin x sin y sin z 
where x, y, z obey the relationship x — y — z,y — x + z and z — y — x and putting the 
condition sin z — 0, the probability expression P^^^^i,^ can be written as 

P,,^,^ = -4 (Pe 



)sin2A;i. (3) 
The condition sin 2; = implies 

I A21 1= invr where n is positive integer. (4) 
Similarly, under the condition sinx = 0, the probability expression Pu^-^v^ can be written as 

P,,^,^ = -4 (Pe 



)sin2A;2 (5) 

and here the condition sin a; = implies 

I Ag^ 1= iriTT where n is positive integer. (6) 

The condition in eq.(4) or in (6) essentially either fixes the length of baseline or the energy 
of the neutrino beam. To determine elements of U' and A^^ we shall use the perturbative 
approach for small and large ^13 seperately. 

We discuss in brief the perturbative approach here. The diagonal neutrino mass matrix 
is approximately given by 

m ^ Arn^idiag(Q^a^l). (7) 
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The effective Hamiltonian induced by interaction of matter with neutrinos is written in weak 
interaction basis as 

H ^ -^R23MRl (8) 

where 

M = i?i3i?i2 R\2R\3 + dtagiA, 0, 0) + 4^ { X* B C \ i?23 ■ (9) 




In equation (9) 



_ 2E^/2GFne _ 2Eeu y _ ^Ee^ ^ _ 2Ee22 ^ _ 2Ee23 ^ _ 2^633 .^q. 

where A is considered due to Standard model interaction of neutrinos with electron. 612 , 
^13, 622, ^23 and 633 are considered due to NSI of neutrinos with matter [9] (e.g, in R violat- 
ing Supersymmetric Models neutrinos may interact with down type quarks through squark 
exchange and may interact with electron through slepton exchange [12] and in Minimal 
Supersymmetric Standard Model with right-handed neutrinos [13] through lepton number 
violating interactions accompanied with neutrinos). We consider magnitude of B,C,D,X,Y 
due to NSI not higher than a due to various experimental constraints [14]. In (9), ( * ) is 
denoted for complex conjugation. In (10), Gp is the Fermi constant and rie is the electron 
number density. 

The mixing matrix U' can be found out as U' — R23 W. Here, W is the normalized 
eigenvectors of AmliM/{2E) calculated through perturbative technique. We follow the 
technique adopted in [15] for Standard Model interactions. Let us consider the case where 
only ei2 and 613 are present as NSI and where sin ^13 is small and of the order of a or less. 
M can be written as M = M(°) -|- M^^^ -|- M^^^ where M* contains terms of the order of a*. 
Then we can write 

M(°) = ^ diag{A, 0, 1), M(^) = ^ ( acj. 



2E •^v ' ' 2E 

A 2 / *13 

Ami ' 



a" 



-e '^aci3si2Si3\ 

M(2) = ^[ -e-^aci2.i2.i3 (11) 

-e'^aCi3s'l2Sl3 -e*^Q!Ci2Sl2Sl3 -s?3 / 



where 



a = C23>^ + e *''ci3Si3 + XS23, b = C23X + Ci2Ci3aSi2 - YS23 ■ (12) 
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The eigenvalues of H upto second order in a are 



/ 2 
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A), 



(13) 



In the same way we can calculate the eigenvalues keeping NSI in 23 block. Using eqs. (4) and 
(13) and putting ei2 and eis to zero one obtains the earlier known magic baseline condition 
[5] in presence of only SMI. For small sin ^13 < a this condition is given by 



L = 2n7r/(\/2G'Fne). 



(14) 



The probability Pi^^-^w^, of oscillation expression after using the basehne condition in (4) for 
such small sin ^13 is (upto order a^) 



'23 



(1-^)^ 



ar sm^— -^^^ 
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(15) 



One can see from this probability expression that \a\ does not contain CP violating phase 
5 when we consider only SMI but it does contain 5 when we keep NSI terms X and Y 
in a. This means that the magic baseline condition (14) is not valid when there is 612 
and as NSI. For brevity, we are not showing the detailed calculation of obtaining the 
probability expression if €22, £23 and 633 is considered. However, one may note that in such 
cases probability P{ve ^n) upto order is same with that in presence of only SMI. Any 
correction due to NSI is present in higher order of a only. So for small sin^^i3, the magic 
baseline condition in (14) is valid when NSI is present only in 23 block of M in (9). 

However, the above conclusions related to magic baseline condition change if we consider 
large sin ^13 > a. Let us discuss the pcrturbative approach for large ^13. This was considered 
in [16] earlier for SMI only. Here we use it for both SMI and NSI particularly in the context 
of magic conditions. We consider NSI in 12 and 13 elements. Then M in (9) can be written 
as M = M(°) + M(i) + M(2) where 
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where 



2E 



ai = C22X - e *'^Ci3Si3Q;Si2 + ^S23- 



(16) 



(17) 



The eigenvalues upto second order correction in a are 
/ 2 
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where 



X= (1+^2-2^008 2^13)^/^ 



(18) 



(19) 



In the same way we can calculate the eigenvalues keeping NSI terms in 23 block also. Using 
(18) the condition (4) may be written for large sin ^13 > a as 



L = 8En7r/ {Am^i (A + l-x)} 



(20) 



Unlike baseline condition in (14) the condition in (20) depends on ^^13 (in x) and energy and 
also this does not give S independent probability as discussed below. For such large sin 613 
the probability P^^^v with baseline condition in (20) is (upto order a^) 



P{ue i^n) ~ -4:Re[Z] sin^ — ^ 2Im[Z] sin 2- 



8E 



8E 



(21) 



where 
Z 



'23 



+ 



;i-a 
C23S23 



a^ek' , aiek' \c2sY + s^sXl'^k' , \c23Y + S2sX\'e 



+ 



X X 
^{C23l3* - S23l*)ik 

{A + l + x) 



+ 



x^ 



x^ 



(22) 



and 



^ = C C23 + -B S23 ; 7 = C23 + C S23 ; C = (-^ + COS 2^13 + a;)cosec2^i3 ; 
k = l/[l + (-A + cos 2^13 + a;)2cosec22^i3]^^^ 



(23) 
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This probability is not independent of 5 due to the presence of Oi in Z. So it is not possible 
to get magic baseline condition (resulting in 5 independent probability P^^^u^ upto order 
a^) for large sin ^13 > a with or without NSI in any elements in M. 

However, if we consider some magic condition on neutrino energy then it is possible to 
get 5 independent probability P^^^u^ for both small and large ^13 and also with and without 
NSI. Using condition (6) and considering P^^^u^ upto order the magic energy condition 
for small sin ^13 < a is written as 

E = AmlJ (±4mr/L + 2\/2GFn^ . (24) 

Using the above energy condition for such small sin ^13 with NSI terms (upto order a^) 

P{^e ^n) ~ -^\c2zX + Ci2Ci3Q;si2 - Fsssl sm^ ( — ^ — j ■ (25) 

With X = y = this corresponds to Standard Model result. Unlike (15) this is independent 
of 5 even with NSI terms. This is one important advantage of using magic energy condition 
instead of magic baseline condition even for sin^^is < a . However, as the condition is on 
energy it might be useful to consider monoenergetic neutrino beam as source [10, 11] to study 
such 5 independent probability. The other alternative way to study such probability might 
be to consider very small energy bins for neutrino energy which satisfies approximately the 
above energy condition in (24). 

Using condition (6) and considering P^^^i,^ upto order a^, the magic energy condition 
for large sin Oi^ > a is written as 

E = Am^i,5 cos 2^i3L^ (^-1 ± ^1 + Q/R^ / (2Q) (26) 

where 

Q = {2m:f - S^L^; R = -S^L^cos^ 2^13; S = v^GpUe. 

Although this condition depends on 613 but with presently allowed values of 613 [2] this 
dependence is not so significant as shown later in Figure 1 in which E satisfying condition 
(24) overlaps on E satisfying condition (26). It is important to note here that unlike magic 
baseline condition (20,) this energy condition (26) results in 5 independent probability as 
shown below. Using condition (26) and considering NSI terms and sin 613 > a, the probabihty 
Pu^^v^ is written as (upto order a^) 

^(^^ ^ " (i+e)(A+i-.)^ (^^^8^ ) ■ ^''^ 
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Unlike the energy condition (26), the probabihty of oscillation depends significantly on ^13. 
This probability is 5 independent with or without NSI. 

Due to present ambiguity in the sign of Amli, with (+) sign to Amg^ for normal and 
with (-) sign for inverted hierarchy of neutrino masses, the energy conditions will be different. 
Apart from hierarchy sign there is further consideration of choosing signs in energy conditions 
as seen in (24) and (26). The requirement of positive energy allows considering both of those 
signs in the energy conditions for normal hierarchy and considering only (-) sign for inverted 
hierarchy. For normal hierarchy for (-) sign in the conditions, L > but there is no 

bound for (+) sign. For inverted hierarchy L < Due to singularity at L = 

neutrino energy E satisfying energy conditions are not possible at magic baseline length. 

Using magic energy condition one could resolve ambiguities in S — 613 (which usually 
happens for non-magic neutrino energy) as the probability is 5 independent. Furthermore, 
using energy condition one may also try to find the neutrino mass hierarchy. For illustration, 
let us consider say nature admits normal hierarchy and sin ^13 < a and for simplicity say 
NSI is absent. Certain neutrino energy has been fixed by energy condition in (24) with 
appropriate choice of n value for which monoenergetic neutrino will be feasible in experiment. 
Under such conditions the probability in (25) is independent of ^13 upto order a^. So the 
probability has fixed value and normal hierarchy could be verified by experiment upto order 

from the number of /i events observed at the detector. This number in general, differs 
from that which one could have obtained for inverted hierarchy in this case. The reason is 
that, for inverted hierarchy the same neutrino energy will not correspond to magic energy 
anymore. In fact, in this case, one can show that the difference in probability of oscillation 
Pve^v^'"'^^ for normal hierarchy with magic energy condition with that for inverted hierarchy 

P,I(non— magic) ..i , . i-j.- • 

i^e->fM Without any magic energy condition is 

pNimag^c) _ plinon-magic) ^ -—1— sin^ 2^13 sin^ ^23 sin^ ( il + ^l^I^l^^ 

+ -rT— — TT« sin 2^13 cos ^13 sin 2^12 sin 2^23 X 

^4(1 + A.) 

( , AmtL\ . f AAm\-.L\ . 
cos (5 sm . „ ^ sm 



V AE ) \ AE ) \ AE J 



which does not vanish in general, for any value of 6 unless sin ^13 vanishes. However, after 
fixing L and E for the experiment one may check whether this difference vanishes or not. 
So for normal hierarchy, in general, one is supposed to get different number of fi events at 
the detector than that for inverted hierarchy and neutrino mass hierarchy may be resolved 
for < sin 613 < a. If the number of fj, events does not match with the expected one for 
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normal hierarchy at magic energy, one may try the magic energy for inverted hierarchy in 
the experiment. Similar to the above expression, one can show that Pj}!^"^*'^^ _ pN(non-magic) 
will not vanish in general, for any value of S unless sin ^13 vanishes. 

For larger sin ^13 > a, the above difference has more complicated form. Same method 
can be adopted in this case also to resolve hierarchy, provided that the difference of P^^^v^ 
(for normal hierarchy) and -Pj^-^.^,^ (for inverted hierarchy) does not vanish for L and E value 
considered in the experiment (where E could be magic energy for either normal or inverted 
hierarchy). In case, it vanishes either for particular combination of sin ^13 or 5 one may 
consider for the same baseline, a different neutrino magic energy by changing n value in the 
magic energy condition for which such difference may not vanish. 

We now, illustrate the use of such magic energy conditions to get 5 independent proba- 
bility and discuss the experimental feasibility. One is required to fine-tune the energy of the 
monoenergetic neutrino beam. In the electron capture facility as discussed at the beginning, 
the neutrino energy may be fixed by appropriately choosing the boost of the ion source. We 
consider the monoenergetic neutrino beam for the ion type with neutrino energy E^ at 

rest given by 1.4 MeV as suggested in ref. [10]. We have chosen Lorentz boost 7 in the range 
of 90 -195 such that the magic energy condition is satisfied. The neutrino energy E is fixed 
in the forward direction by the boost as £^ = £'^7. We have assumed flux of 10^^ neutrinos 
per year. We are considering a baseline of length 650 Km from CERN to megaton water 
Cerenkov detector possibly located at Canfranc in Spain. For such basehne the constant 
matter density has been approximated to be 4.21 gm/cc. For our subsequent sensitivity 
analysis of oscillation parameters we mention here the detector characteristics also for such 
experimental set-up. [7]: (a) Fiducial mass = 500 Kton (b) Detection efficiency (e ) = 50 % 
(c) Charge identification efficiency (/g) = 95 %. 

In Figure- 1 we have shown the energy versus length of baseline satisfying magic energy 
conditions. Condition (24) has been considered for any sin ^13 < a and condition in (26) for 
^13 = 5, 8, and 12 degrees. However, it is seen from the figure that the plots with different 
^13 are almost overlapping with each other indicating very small change in energy E with L 
due to variations of unknown parameter ^13. In plotting Figure 1, instead of n = 1 we have 
considered n = 2 in the energy conditions (24) and (26) so that for the above-mentioned 
baseline of length 650 Km, the magic neutrino energy lies in experimentally feasible range. 

Finally we discuss the sensitivity in measuring the unknown oscillation parameters like 
^13 and Amg]^ in the experimental set-up with monoenergetic neutrino beam [10, 11]. The 
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n=2 



Normal Heirarchy 
Inverted Heirarchy 



Figure 1: Plots for neutrino energy (E) versus length (L) of baseline using energy conditions (24) 
for any sinOi^ < a and (26) for ^13 = 5, 8, 12 for n = 2 for both hierarchies. 



number of n events expected at the detector due to ^'/^ oscillation is given by 

N^ = TnnIeeE <f){E) a,^ {E)P,^^,^ {E) (28) 

where T = time period , rin = number of target nucleons, 4>{E) = flux, a^^^E) = detection 
cross-section. As we are considering monoenergetic neutrino beam so the energy resolution 
function on which normally number of events depends, may be considered to be effectively 
1 and as such is not mentioned in eq. (28). In figure 2 we have shown for both hierarchies 
the variation of the number of events expected for a priod of 5 years with ^13 in absence 
of any NSI and with NSI respectively. The number of fi events are quite large as neutrinos 
have some fixed energy instead of Gaussian distribution of energy. Number of events with 
inverted hierarchy is found to be slightly higher than that for normal hierarchy for different 
values of ^13 with or without NSI. 

In plotting figures 4 and 5 we define Xtotai 

2 _ /AT^^pi _ Ar*^\ 2 /|Am2j - |Am|i(trae)|y /sin^ 2^23 - sin^ 2^23 (trwe) V 
" V AiV ) + I ^A^) ) + 1^ a(sin^2^^23) J ^ ^ 

where N'^^^^ and A^*'^ stands for experimental and theoretical number of /i events respec- 
tively at the detectors and the error in A^*^ = AA^ = AApert + AA"q,2. Here, ANp^rt and 
AN^-i are the differential change in A^*^ considering perturbative expression of Pi^^^y^ in (27) 
and the perturbative error of order in (27) respectively. The AA"q,2 takes care of maxi- 
mum possible correlated error due to any true value of 5 and we have assumed democratic 



10 



29 




0.05 0.10 0.15 0.20 

01 3 (radian) 

Figure 2: Number of ji events versus mixing angle ^13 with 5 years running without NSI (i.e 
X = Y = 0). 



form of correlation matrix. Unlike magic baseline condition for magic energy condition the 
perturbative expression of probability depends on cos^ ^13 instead of sin^ 26^13 and thus re- 
sulting in large number of events (~ 10^) even for small di^. So unlike [7] instead of Poisson 
distribution, we have considered Gaussian distribution for error in A^. In evaluating 2nd 
and 3rd terms on right side of eq.(29) we have considered Amli{true) = 2.5 x 10~^ eV'^, 
sin^ 26^23 (trwe) = 1.0, cr(Am|]^) = 1.5%, a (sin^ 26'23) = 1%. Based on various experimental 
data set [7] the following 3a constraints on the following oscillation parameters have been 
considered : 

2 X 10"^ eV^ < lAm^J < 3.2 x 10"^ eV^; 
35.67° < 023 < 55.55°. 

In figure 4 and 5 we have shown the precision in measuring cos^ 613 which appears in the 
expression of Py^^^^ satisfying magic energy condition. This is defined by 

„ cos^ ^^13 (mm) — cos^ 6'i3 (max) ^ , , 

Precision = — , . f 100% (30) 

cos"^ Ui3\min) + cos"^ Ui2,\max) 

in which cos^6'i3 (mm) and cos^ 6'i3 (max) are the smallest and largest values respectively of 
cos^ ^13 at the given confidence level. We have shown the precision at la and 3(T level with 
and without NSI for both the hierarchies as mentioned in the figures. In finding precision 
in measurement of cos^ ^13 in presence of NSI, we have assumed that the strength of NSI 
couplings are known from some other experiments. The number of muon events at the 
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Figure 3: Number of n events versus mixing angle ^13 with 5 years running with NSI X and Y. 
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(a) (b) 
Figure 4: Precision in the measurement of cos^ ^13 without NSI expected with 5 years running at 

la in (a) and at 3a in (b). 



detector changes when such NSI are included and thus changes the level of precision in 
measurement of cos^0i3. The precision as defined above is better for its lower values. In 
general in both figures 4 & 5 the precision initially deteriorates with increase of cos^ 6'i3 
values and after reaching some limiting values of cos^ ^13 the precision starts improving and 
finally for larger values of cos^ ^13 again it deteriorates. However, in some cases the change in 
precision is insignificant with the variations of cos^ ^13 as seen in figures 4(b) and 5 (c) & (d). 
From figures the differences in precision for normal and inverted hierarchies are particularly 
found at la level except in figure 5 (c) where non-zero value of Y as NSI has been considered. 
Comparing three cases - (a) no NSI (b) NSI with only non-zero X (c) NSI with only non-zero 
Y, it is seen from figures 4 & 5 that the precision in measurement of cos^ ^13 is best for (b) 
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Figure 5: Precision in the measurement of cos ^13 with NSI expected with 5 years running at la 
in (a) and (c) and at 3a in (b) and (d). 



and the precision is better for (a) than that for (c). Actually more the number of /i events 
the better is the precision as these cases may be seen from figures 2 & 3. 

As concluding remarks we mention that to get P^^^^^ almost independent of unknown 
CP violating phase 6 one may consider either magic baseline condition on the length of 
baseline or the magic neutrino energy condition. However, there are some disadvntages in 
considering the magic baseline condition which are not present when magic energy condition 
is considered. The magic baseline condition exists only for small sin6'i3 < a. Also this 
condition exists when NSI are considered in only 23 block of effective neutrino mass matrix 
M (as Pu^^u^ upto order is independent of NSI in 23 block). However, magic baseline 
condition is not possible if NSI terms are present in 12 , 13 elements of M. For large 
sin 6*13 > a using magic baseline condition in (20) it is not possible to get 6 independent 
probability P^^^y^ with or without NSI. Magic baseline condition will also depend on neutrino 
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energy for sin ^13 > a. Also to place neutrino detector at a location satisfying magic baseline 
condition may not be always feasible. Some of the drawbacks mentioned above in considering 
magic baseline condition may be overcome by considering the condition on neutrino energy. 
Under magic neutrino energy condition for both small and large sin ^13 and also with or 
without NSI one gets 5 independent probability Pn^^^^. 

Using energy condition there is scope to find out the hierarchy of neutrino masses and 
to obtain overall good precision in the measurement of cos^ ^13 over the allowed range of 
^13 as discussed earher. Depending on the presence or absence of NSI, the number of ji 
events could significantly differ as shown in Fig. 2 & 3 and could signal the presence of 
new physics. If the experimental data indicates the presence of NSI then to find NSI as 
well as ^13 one may consider changing the magic neutrino energy by changing the Lorentz 
boost in the same experimental set-up. Then NSI terms like 612 and 613 as well as ^13 may 
be known from P^^^^i,^ in (25) for sin ^13 < a or (27) for sin ^13 > a after matching those 
probabilities with experimntal data on number of // events using eq. (28). In long baseline 
experiments, monoenergetic neutrino beam as source with neutrino energy satisfying magic 
energy condition, could be highly useful in future advanced precision measuments of neutrino 
oscillation parameters, in resolving hierarchies of neutrino masses and also in searching NSI 
of neutrinos with matter. 
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